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Abstract 

Experiments with H~ atoms are proposed in order to study the nuclear 
interaction of S hyperons. The production of H~ in the (K~ , K + ) reaction, 
the H~ stopping in matter, and its atomic cascade are incorporated within 
a realistic evaluation of the results expected for H~ X-ray spectra across the 
periodic table, using an assumed H-nucleus optical potential V op t- Several 
optimal targets for measuring the strong-interaction shift and width of the X- 
ray transition to the 'last' atomic level observed are singled out: F, CI, I, Pb. 
The sensitivity of these observables to the parameters of V op t is considered. 
The relevance of such experiments is discussed in the context of strangeness 
—2 nuclear physics and multistrange nuclear matter. Finally, with particular 
reference to searches for the H dibaryon, the properties of E~(f atoms are also 
discussed. The role of Stark mixing, its effect on S and P state capture of 
H~ by the deuteron, together with estimates of the resulting probability for 
producing the H dibaryon are considered in detail. 
PACS: 15.20.Jn, 36.10.Gv 
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I. INTRODUCTION AND BACKGROUND 



Very little is established experimentally or phenomenologically on the interaction of H 
hyperons with nuclei. Dover and Gal analyzing old emulsion data which had been 
interpreted as due to H~ hypernuclei, obtained an attractive S-nucleus interaction with a 
nuclear potential well depth of Vq = 21 — 24 MeV. This range of values agrees well with 
the theoretical prediction by the same authors || for H in nuclear matter, using model D 
of the Nijmegen group H to describe baryon-baryon interactions in an SU(3) picture, in 
contrast with the S- nucleus repulsion obtained || using model F [|J. Similar predictions 
were subsequently made with more detailed G matrix evaluations by Yamamoto et al. |5]|| 
who argued for a considerable A dependence of Vq"\ such that the well depth for light and 
medium weight nuclei is significantly lower than for heavy nuclei where it approaches the 
value calculated for nuclear matter. It should be noted, however, that the predictions of the 
Nijmegen model D for Vq are extremely sensitive to the value assumed for the hard-core 
radius. Nevertheless, the confidence in the predictive power of model D for this sector of 
strangeness —2 hypernuclear physics, at least qualitatively, is to a large extent due to its 
success in yielding the attractive AA interaction necessary to reproduce the (so far) three 
known AA binding energies (see Ref. || for a review of these calculations). 

If the interaction of H hyperons with nuclei is sufficiently attractive to cause binding, 
as has been repeatedly argued since the original work of Dover and Gal |lj, then a rich 
source of spectroscopic information becomes available and the properties of the in-medium 
EN interaction can be extracted. Bound states of 5 hypernuclei would also be useful as a 
gateway to form double A hypernuclei f7],§]. Finally, a minimum strength for V of about 
15 MeV is required to realize the exciting possibility of strange hadronic matter ||; where 
protons, neutrons, As and Ss are held together to form a system which is stable against 
strong-interaction decay. 

Some new information on the S~ nucleus interaction has been recently reported from 
(K~, K + ) counter experiments at the KEK proton synchrotron. Fukuda et al. |Hj have 
shown fits to the very low energy part (including the bound state region) of the S~ hy- 
pernuclear spectrum in the l2 C(K~, K + )X reaction on a scintillating fiber active target 
(experiment E224), resulting in an estimate of Vq^ between 15 to 20 MeV. The experimen- 
tal energy resolution of about 10 MeV in this experiment was too poor to allow identification 
of any bound state peak structure which could have given more definitive information on 
the well depth. A somewhat cleaner and better resolved spectrum has been recently shown 
1 1] from the Brookhaven AGS experiment E885, but no analysis of these data has yet been 



reported. An earlier KEK experiment (E176) gave evidence for three events of stopped S~ 
in light emulsion nuclei, each showing a decay into a pair of single A hypernuclei. The 
first two events HT2| , |l"3f are consistent energetically with a S~ atomic state in 12 C bound by 
_B=-( 12 C) = 0.58 ± 0.14 MeV. However, this value could only be ascribed to capture from 
the IS state which is estimated to occur in less than 1% of the total number of captures. 
This binding energy is distinctly larger than the calculated value B~-( 12 C) < 0.32 MeV for 
the 2P state, for a wide range of strong-interaction potentials. Moreover, the S~ capture 
probability in 12 C from P states is a few percent at most. The most likely capture in 12 C, 
as discussed in Sec. I ITB, occurs from atomic D states. The calculated binding energies 
of the atomic 3D states for C, N, O emulsion nuclei are given in Table III where it is seen 



1 



that binding is essentially by the Coulomb potential. The two examples for binding in the 
presence of a strong H-nucleus potential are for the tp potential used in Sec. || and |TT| with 
the parameter b = 0.25 + i0.04 fm (potential 1) and b = 0.19 + i0.04 fm (potential 2), 
corresponding to Vq =20.5 and 15.6 MeV, respectively, in 12 C. (The value used for Im5 
is discussed in Sec. |III A| and fV\j . It is seen that the sensitivity of the binding energies 
and strong interaction widths to the S-nucleus strong interaction is, for these examples, 
of the order of 100 eV, substantially smaller than a typical error of 100 keV incurred in 
emulsion work. We point out that there exist alternative interpretations of these two events 
as captures on 14 N, with binding energies consistent with the calculated value listed in the 



table, for example [EJ, 5=~( 14 N) = 0.35 ± 0.20 MeV. Furthermore, a likely interpretation 
of the third event |3 is due to capture on 16 0, with B E -( w O) = 0.31 ±0.23 MeV. Clearly, 



whereas these emulsion events are consistent with capture from 3D atomic states, they are 
useless as a source of information regarding the S-nucleus interaction. 

One clearly needs an alternative source of information on the H-nucleus strong inter- 
action. Such an alternative source is the measurement of X-ray energies from transitions 
between low lying levels of H~ hadronic atoms. The experimental accuracies of the proposed 
measurements are such that meaningful information on the H-nucleus interaction is likely 
to be obtained. In arguing the case for doing experiments with H~ atoms, we follow the 
example of atoms which, as was recently shown [15], gives rise to meaningful information 
about the £ nucleus interaction, particularly in the absence of any systematic evidence for 
bound E hypernuclei. 

Experiments with stopped H~ hyperons had been proposed by Zhu et al. |16| and by 
Akaishi and collaborators |l7j in order to produce some of the lightest AA hypernuclei, 
Ai YHe and AA H (if the latter is particle stable) respectively, by looking for a peak in the 
outgoing neutron spectrum in the two-body reaction 
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These proposals motivated the AGS experiment E885 |nj on 12 C, using a diamond target 



to stop the 5~ hyperons resulting from the quasi-free peak of the (K~, K + ) initial reaction. 
Finally, stopping S~ hyperons in deuterium has been used in the AGS experiment E813 to 
search for the doubly strange H dibaryon through the reaction (E~d) atom — > Hn, as reviewed 
recently by Chrien |18| 



In Sec. m we consider the experimental features of measuring X-ray transitions in 2~ 
atoms, and in Sec. [I il] we discuss possible target nuclei, using a tp optical potential. Al- 
ternative choices, using density dependent optical potentials, are discussed in Sec. [TV} The 
special case of H~ atoms of deuterium, which is connected to recent searches for the H 
dibaryon, is discussed in Sec. M. Section summarizes the present work. 



II. EXPERIMENTAL CONSIDERATIONS 

In considering the feasibility of measuring strong interaction effects in S~ atoms, we are 
guided by the successful observation and measurement of strong interaction effects []19|-|2ll 
in E~ atoms. The S~ and the £~ hyperons have very similar masses and lifetimes, namely 
1321.32 and 1197.34 MeV for S" and E~ respectively, and 0.1642 and 0.1482 nsec for S~ and 
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£~ respectively. Consequently, major differences between the experimental X-ray counting 
rates should result mostly from differences between the production mechanism of the two 
kinds of hadronic atoms. In the E~ case the production mechanism is the K~p — > £~7r + 
reaction at rest. The outgoing S~ has a kinetic energy of 12.4 MeV which means that it 
stops in a heavy target such as Pb or W in less than 1CT 11 sec and the outgoing pion has 
a kinetic energy of 83.2 MeV which makes it quite easy to detect. For the H~ particle the 
situation is less favourable. The production mechanism is the K~p — > E~K + reaction 
which has a total cross section larger than 100 /A) only in the kinetic energy range for K~ 
of 1.0 to 1.8 GeV (see results summarized in Ref. |]T|). Consequently the majority of the 
produced S~ will be in the energy range of 100 to 300 MeV, leading to stopping times in 
heavy targets of 0.01 to 0.3 nsec respectively. It is therefore essential that the S~ particles 
slow down in a heavy degrader and not in the liquid hydrogen used for the production where 
the stopping times are much longer. A laminar target structure similar to that used in £~ 
experiments seems the most suitable. Without detailed calculations for specific target 

and degrader configurations it is only possible to estimate that up to 50% of the produced 
2~ may survive the slowing down process. Detection in coincidence with the outgoing K + 
is a must in such an experiment. The energy of this K + will be about 500 MeV. 

When interested in the measurement of strong interaction effects in hadronic atoms it 
is also necessary to consider losses during the atomic cascade process. Here, again, it is 
possible to make use of the fact that S~ atoms have been successfully observed, and that 
the S~ and the S~ hyperons have very similar masses and lifetimes. We have performed 
calculations of the atomic cascade for S~ atoms of Pb, in order to compare the predicted 
absolute X-ray yields with experiment JJl[ . 



The cascade program used here was based on that originally written by Hiifner |23j for 
muonic atoms and later modified to include strong interaction effects. The hadronic atom 
is assumed to be formed in a state of large radial quantum number n. To accommodate 
departures from a purely statistical distribution the initial population used is of the form 

p(l) = (21 + l)e al (2) 

where for a = the statistical distribution is obtained. The atom initially de-excites by 
Auger transitions and later by radiative transitions with the emission of X-rays. Finally in 
a state of low angular momentum I the hadron is absorbed by the nucleus and the X-ray 
cascade terminates. The probabilities for Auger and radiative transitions are calculated 
with the usual electromagnetic expressions ||24|| . For the strong interaction the input to the 



calculations are the measured or calculated absorption widths of low n circular (n, I = n — 1) 
atomic states which are then scaled [[ZIJ to obtain widths for non-circular states of a given 
I and increasing n. The probability that the hadron may decay during the atomic cascade 
is also included in the calculation. 

For the strong interaction widths we used values predicted by the density dependent E 
nucleus potential which reproduces all the available data ||15|| . From a fit to the measured 
relative X-ray yields |2]J we obtain the value a = —0.056 ± 0.020 which corresponds to a 
small departure from a purely statistical distribution. The predicted absolute yields for the 
three transitions of interest are shown in Fig. [l| where it is seen that their variation with 
a over its fitted range of uncertainty is not very large. Note that the 11M — > 10L and 



10L — > 9K transitions in £ Pb have been observed [21]. Next we calculated the yields for 
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2~ atoms of Pb for a range of values of a. Strong interaction widths of the relevant levels 
were calculated with a tp potential with 6 = 0.25 + i0.04 fm as discussed below. Figure |2| 
shows the calculated absolute yields, which turn out to be very similar to the corresponding 
values for the E~ Pb atom. In both cases the 13 — ► 11 transition is also shown because its 
energy is very close to that of the 10 — > 9 transition and the two will have to be separated. It 
is therefore concluded that S~ atoms could be observed if the production yield and detection 
efficiency are not too small compared to the corresponding values for S~ atoms. 



III. SELECTION OF TARGETS 
A. Potentials and criteria 

When selecting targets for possible experiments on S~ atoms, it must be assumed that 
such experiments will probably not be feasible on more than very few targets, and one must 
therefore ask whether it is at all likely that useful information on the interaction of with 
nuclei will be obtained from the resulting rather limited range of data. Useful hadronic-atom 
data normally consist of the strong- interaction shift and width for the 'last' level observed 
plus, occasionally, the relative yield for the 'upper' level, as discussed below. It has been 
shown very recently |^B| that the main features of the interaction of K~ and £~ with nuclei, 
as found from analyses of all the available data, may in fact, be obtained by analyzing a 
small fraction of the available hadronic atom data, if the target nuclei are carefully selected. 
A key point here is to have target nuclei over as wide a range of the periodic table as possible. 
This observation suggests that experiments on H~ atoms may provide useful information. 

In order to have some idea of the expected strong interaction effects in H~ atoms and 
on the expected yields of X-ray transitions, it is necessary to adopt some form of an optical 
potential which will describe, at least approximately, the interaction of the S~ hyperon with 
nuclei and its dependence on the nuclear mass. The so-called tp potential |27j in its simplest 
form is given by 

2pV opt (r) = -4ir(l + j^)b p(r) (3) 

where p is the S-nucleus reduced mass, M is the mass of the nucleon and p(r) is the nuclear 
density, normalized to the mass number A. The parameter bo is a complex parameter, and 



is usually related to the hadron- nucleon scattering length |[27|| . More refined potentials have 
been used with other hadronic atoms and will be discussed in the following section. For 
the present purpose we select the value of Re6 = 0.25 fm which yields a potential depth 
of about 20 MeV inside nuclei, and Im6 = 0.04 fm, yielding for the imaginary potential a 
depth of about 3 MeV. Whereas the real potential may be regarded as 'typical', according 
to the discussion in Sec. |, the imaginary potential is about twice as large as estimated || 
in model D. Reducing the imaginary potential will only cause the calculated widths of the 
states to decrease by roughly the same proportion, and the relative yields (see below) of 
transitions to become larger. This will not, however, change the last observed atomic level. 



A further comment on the choice of the imaginary potential is made in Sec. |TV . 



In choosing criteria for the suitability of a transition as a source of information on the 



H nucleus interaction, we are guided by experience with other hadronic atoms [27,28] and 
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select X-ray transitions (n + 1, 1 + 1) — > (n, i) between circular atomic states (n = / + 1) 
with energies greater than 100 keV, where the strong interaction shift for the 'last' (n, I) 
level is at least 0.5 keV and the width less than about 10 keV. The 'upper' level relative 
yield, defined as the ratio of the intensity of the (n + 1, 1 + 1) — > (n, I) X-ray transition to 
the summed intensity of all X-ray transitions feeding the (n + 1, 1 + 1) state, is also required 
to be at least 10%. The 'upper' level relative yield p8| is given by the ratio of the radiative 
width for the (n + 1, 1 + 1) — > (n, /) transition to the total width of the (n + 1, 1 + 1) level 
feeding the (n, /) level. A measurement of the relative yield enables the strong interaction 
width of the 'upper' (n+ 1, 1 + 1) level to be determined, in addition to deducing the strong- 
interaction level shift and width for the 'last' (n, I) level from the X-ray transition energy 
discussed above. Altogether, these are the strong-interaction data provided by measuring 
atomic X-ray spectra. 



B. Results 



Strong interaction shifts and widths of S~ atomic levels have been calculated using the 
above optical potential for a large number of nuclei. As the overlap of atomic wavefunctions 
with nuclei vary smoothly with charge number, it is to be expected that generally shifts, 
widths and yields will vary smoothly along the periodic table. Figure |3] shows calculated 
widths and 'upper' level relative yields for the 9K level in S~ atoms for several heavy targets. 
Calculated shifts (not shown) for the chosen potential are generally equal to or up to 30% 
larger than the widths, as long as both are smaller than 10 keV. It is seen from the figure 
that Pb is a suitable target for the experiment. Figure |] shows similar results for the 71 
state in medium-heavy S~ atoms and it is seen that a suitable target may be found near Sn 
or I. The dashed lines in this figure are obtained by reversing the sign of the real potential 
used for calculating the solid curves. It is seen that in such a case the range of suitable 
targets will move to between I and Ba, where the strong interaction width and relative yield 
are more acceptable. The sign of the strong interaction shift will be reversed in this case, 
but it has no experimental consequences. This exemplifies a general property of hadronic 
atoms, which are dominated by the Coulomb interaction, namely, that large variations in 
the strong interaction potential will move the proposed targets only a few units of charge 
along the periodic table. 

Figure |5| shows results for the AF state of S~ atoms, where it is seen that for a Si target 
the effects could be too small to measure whereas for Ca the width could be too large and the 
relative yield too small. In this region a CI target may be appropriate, perhaps in the form 
of the liquid CCI4. More detailed results for CI are shown in Fig. ^| where the sensitivities 
to assumptions regarding the optical potential are also typical of results for other targets. 
The solid curves connect points obtained within the tp potential Eq. (Q) for fixed values of 
Re6 , listed above the lines. The four points along each line correspond to values of Im5 
from 0.05 fm down to 0.02 fm in steps of 0.01 fm. Departures from this tp potential are 
represented by the dotted lines, calculated from phenomenological density dependent (DD) 
real potentials similar to those found from analyses of experimental results for S~ and K~ 
atoms ||27fl . The imaginary part of the potential is of the tp type and the points along the 
dotted lines correspond to the same values of Imfe as above. The real potentials in these 
calculations are similar to the real potential for S _ atoms |27|| , having an attractive pocket 
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about 5-10 MeV deep outside the nuclear surface, with a repulsive potential of about 20-30 
MeV in the nuclear interior. The results in the figure serve only to illustrate the expected 
range of strong interaction effects. If the actual values of shift and width turn out to be 
within the area covered by the lines, these effects will most likely be measurable. 

Figure ^ shows calculated results for 3D states in several H~ atoms, and it is seen that 
F may be a suitable target, possibly in the form of teflon (CF2). To see what could be the 
effects of the carbon in the target we have also performed full atomic cascade calculations 
for C and F, and the expected X-ray spectrum for teflon is shown in Fig. §. It seems that 
the presence of carbon in the target should not affect the possibility of observing transitions 
in H~ F atoms. 

Results for the tp potential with bo = 0.25 + i0.04 fm are summarized in Table f|, 
which includes several suitable targets across the periodic table where strong interaction 
effects are likely to be measurable. A range of targets is necessary if unusual features of the 
interaction, such as density dependence beyond the tp prescription, are to be unraveled by 
the experiment. 

IV. 3 NUCLEUS DENSITY DEPENDENT POTENTIALS 

In the previous section we have discussed H~ X-ray spectra across the periodic table, 
assuming a H-nucleus strong interaction potential of the form = tp, Eq. ([5]). Specific 
predictions were given, using the value b = 0.25-H0.04 fm, and the sensitivity to departures 
from this value was displayed in Fig. ^ and |6|. Conversely, one may explore the extent 
to which measuring, say, four X-ray spectra (e.g. Table [n]) will determine the H-nucleus 
potential. In the present section we address this question by choosing a one-boson-exchange 
(OBE) motivated H-nucleus potential, use it to evaluate the strong-interaction shifts and 
widths for several targets across the periodic table, and then solve the inverse problem, 
namely determine the best-fit parameters for H nucleus potentials of different forms. 

We start with the HA DD G-matrix interaction YNG derived by Yamamoto et al. || 
from the Nijmegen OBE potential model D |3|: 

3 

v-en (r; k F ) = ( a i + hk F + Cik F ) ex P (~ r V A 2 ) > ( 4 ) 

i=l 

where the nuclear Fermi momentum hp is expressed in terms of the local nuclear density 

1 /3 

p, kp = (3ir 2 p/2) . To simplify matters, while retaining a substantial part of the density 
dependence of v^n, we kept only the a and b terms. These G- matrix parameters are listed in 
Table for the choice r c = 0.46 fm, which belongs to the range of values considered in Ref. 
0. We note that the G matrix depends strongly on the hard-core radius r c chosen in the 
HA sector of model D. The resulting Vsn is repulsive at short HA distances and attractive 
at distances r > 0.8 fm. The H nucleus DD real potential, obtained by folding v^.n with the 
nuclear density p: 

(r; p) = J v EN (|r - r'| ; kp (p)) p (r') d 3 r' , (5) 

is shown in Fig. ^| (G, solid line) for the nucleus Sn. Its depth is about 18 MeV (attractive) in 
medium- weight and heavy nuclei and about half of that in light nuclei, so it falls conveniently 
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within the range of values considered as 'realistic' in Sec. |. However, as is clear by comparing 
the shape of this V^"' (potential G) with the (fitted, see below) tp potential, plotted as a dot- 
dash line in the figure, potential G extends considerably further out than does the nuclear 
density distribution p . The difference in the r.m.s. radius values amounts to about 1.3 fm. 
The main contribution to this difference does not arise just from folding the finite-size v^n 
into p, but is due to the summed effect of the DD terms l V which produce repulsion on 
top of the attractive contribution due to the sum of the 'a' terms. Retaining the 'c' terms 
in Eq. (|j) would have reduced the r.m.s. radius of potential G, but in addition it would 
have increased the nuclear-matter depth to about 30 MeV, a value which currently []TU[ is 
considered to be a gross overestimate. 

For the imaginary part of V^"' we retained, for simplicitly, the tp form with Im&o = 0.04 
fm as before. We checked that replacing it by the corresponding DD form ||, as per Eqs. 
(DID, has a marginal effect on our results. The reason is that since the real part of V^~' is 
substantially stronger than the imaginary part, the main effect of the density dependence is 
due to the real part of . 

The calculated strong-interaction shifts and widths, using potential G, are shown in 
Table |V| for several target nuclei, including the five targets listed in Table |H] which were 
motivated by the tp potential with a similar depth. Where comparison can be made, these 
shifts and widths are larger by an order of magnitude than those in Table [H], due to the 
greater extension of potential G which pulls in the atomic wavefunctions by up to 2 fm 
(judging by the inward shift of the position of the maximum of the atomic wavefunctions 
at about 20 fm). Consequently, with most of the calculated shifts and widths exceeding 10 
keV, the targets (F, CI, I, Pb) will no longer be the optimal ones for experiments with 2~ 
atoms. As discussed in Sec. |I I [B|, one can identify neighbouring targets by going several 
charge units down in each sequence to reestablish optimum conditions. Indeed, according 
to Table |TV|, an optimal set of targets could consist of (O, S, Sn, W). 

For testing the inversion procedure under the extreme conditions set up by using potential 
G, the shifts and widths calculated for these four S~ atoms were considered as pseudo data. 
Standard measurement errors of 20% were ascribed, and randomized values of these pseudo 
data quantities were then derived and used for fitting a given form of the real potential 
(the imaginary part was always taken as tp). Two forms were used: 

(i) The tp form, but with bo determined from the fit. The best fit parameter bo was found 
to be 



Re&o = 0.402 ± 0.022 fm, Im6 = 0.053 ± 0.021 fm, 

with x 2 /F) the \ 2 of the fit per degree of freedom, assuming a relatively high value of 
X 2 /F = 8.8. The real part of this potential is shown in Fig. |^. Note that it is twice as deep 
in the nuclear interior as potential G, owing to the analytic continuation from the surface 
region where the fitted potential must assume larger values than before. However, with such 
a distinct difference between the geometries of the two potentials, the quality of fit produced 
by the tp best-fit potential to the pseudo data due to potential G is rather poor. 

(ii) A phenomenological DD potential of the form Eq. (|3|), where the parameter 6 now 
depends on density: 

bo(p) = bo + B [p(r)/ P (0)r. (6) 
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For the imaginary part, the tp form Eq. (0) was maintained, i.e. ImBo = 0. If the exponent a 
is taken as a = 1/3, in order to simulate the density dependence of the underlying G-matrix 
vzn (Eq. ((D), then the best-fit parameters are 

Re&o = 1-27 ± 0.21 fm, Im6 = 0.019 ± 0.009 fm, 
ReB = -1.05 ± 0.26 fm, x 2 /F = 4.65. 

This potential is similar in depth to potential G, its r.m.s. radius is substantially larger 
than that of the tp potential, and as p —>■ it is about 3 times stronger than the best-fit tp 
potential discussed above. The quality of this fit is quite reasonable. 

If a is allowed to vary, it is found that the value a = —0.7 provides a very good fit to 
the pseudo data: 

Re&o = 0.268 ± 0.033 fm, Im6 = 0.040 ± 0.009 fm, 
ReB = 0.056 ± 0.011 fm, j?/F = 1.77. 

This best-fit phenomenological DD potential is also plotted in Fig. |9|. It is deeper within the 
nucleus than potential G, becoming less attractive than the latter in the immediate surface 
region outside the nucleus, and then (since a is negative) becoming more attractive in the 
more extreme surface region (8 fm and beyond, for Sn). 

It is clear that neither the oversimplified tp potential, nor the phenomenological DD po- 
tentials defined by Eq. (§), are capable of reproducing the oscillatory behaviour of potential 
G in the nuclear surface region. The best-fit phenomenological DD potential probably gives 
a reliable indication of the overall attractive nature of the underlying potential both 
yielding similar strengths of about 10 MeV in the nuclear surface region. It is noteworthy 
that the fitted imaginary part is practically the same for both potentials, probably due to 
the same tp form assumed in the two cases. 



V. H~ DEUTERIUM ATOMS 

In the previous sections we have discussed the possible observation of X-rays from H~ 
atoms formed in nuclear targets and the possibility of obtaining information about strong 
interaction effects in atoms. In the present section we discuss the properties of 
atoms and in particular the formation of the doubly strange H dibaryon through the reaction 
(E~d) atom — > Hn. An experiment to search for this reaction is in progress [p9] . The formation 



of the H particle is detected by observation of the monoenergetic neutron, whilst the S~ 
is tagged by the observation of the K + from the production reaction K~p — > E~K + as 
discussed earlier in Sec. |I[ The probability for H dibaryon formation from "Erd atoms has 
been calculated by Aerts and Dover |30| for reactions involving S and P states of the E~d 



atom as a function of the mass of the H particle. (See Table II and Fig. 4 of Ref. [fjOfl ). 
Here we are specifically concerned with calculating the probability for the H~ to interact 
with the deuteron from an atomic S or P state in the S~d atom. 

For the E~d atom, the atomic cascade proceeds generally as described in Sec. [TT] but 
with an additional complication due to the presence of Stark mixing ||31|| . This gives rise to 
transitions of the type (n, I ± 1) <-» (n, I) which increases the probability of interactions 
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from atomic S states at large values of n. Since the Stark mixing is proportional to the 
collision rate, and hence target density, the effects are largest in liquid hydrogen. As the 
interactions take place at high n, the lifetime of the E~d atom is reduced by the Stark mixing 
and the probability for the S~ to decay during the atomic cascade decreases. 

Some of the early qualitative discussions of Stark mixing [31] for K~p atoms were placed 
on a more quantitative basis by Leon and Bethe [32|. Their calculations were later mod- 
ified and extended by Borie and Leon [33| to a range of exotic-hydrogen atoms such as 
fi~p, ir~p, K~p and pp. Of interest for the present work, calculations using the Borie and 
Leon model have been used to fit measured X-ray yields for K- and L-series X-rays from 



antiprotonic hydrogen (pp) atoms P^fl . In using this model it has generally been found that, 
because of approximations in the calculation of the absolute value of the Stark mixing rate, 
an overall normalisation factor ksrx is required to get a good fit to the experimental X-ray 
data. Typical values of k STK required to fit data for pp atoms are in the range 1.0 to 2.0. 
In the present work calculations are made for ksrx over the range 0.5 to 5.0. 

The present calculations were carried out using the Borie and Leon model [33[] with the 
program used in for calculations on antiprotonic hydrogen atoms, modified for E~d 
atoms and including decays of the H~. The principal input values are the Stark mixing 
normalization factor k$TK and the strong interaction widths Tig and T 2 p for the IS and 2P 
levels respectively. The strong interaction width of the 3D and higher (n,l = n — l)-states 
were all set to zero. The kinetic energy of the S~d atom was taken to be 1 eV |33fl . 

The total strong interaction widths calculated by Aerts and Dover 0] lie in the range 
Tis = 1 to 50 keV for S states and T 2 p = 0.2 to 0.8 eV for P states. In Fig. [HJthe fraction 
of 5* and P state reactions as a function of the Stark mixing parameter ksrx is shown for 
three sets of strong interaction parameters, T 1S = 10 keV, Y 2 p = 0.2 eV (labelled 1), r is 
= 10 keV, r 2P = 0.8 eV (labelled 2) and T 1S = 1 keV, r 2P = 0.2 eV (labelled 3). The 
fraction of S~d atoms which decay is also plotted for case 1; cases 2 and 3 give almost 
indistinguishable values. The results are seen to exhibit only small sensitivity to the choice 
of strong interaction parameters, but considerable sensitivity to the value of ksrx- As the 
amount of Stark mixing increases, the fraction of 5* state capture increases, P state capture 
decreases and the number of H~ decays from the S~d atom is reduced as expected. 

A particular feature of the Borie and Leon model is that it seems to over-estimate the 
amount of P state capture. This problem has been discussed by Batty |J5| for pp atoms. 
He finds that the introduction of an additional normalization factor Kq, used for Stark 
transitions between S and P states only, i.e. nS — > nP and nP — > nS transitions, gives 
significantly improved fits to the X-ray yield data with K = 7.6 and much reduced values 
for the fraction of P state annihilation in liquid hydrogen. 

Here we use the same model and again fix Kq = 7.6. There is no particular reason 
for this choice of value for K except that it gives a good fit to the pp atom results. The 
resulting values for the fraction of S and P state capture and of E~ decays are plotted in 
Fig. [TT] for the same values of the strong interaction parameters as used in Fig. [L(] which is 
equivalent to the case with K = 1.0. Setting K = 7.6 as in Fig. decreases the P state 
fraction and increases the 5* state fraction, as is to be expected for the increased amount of 
Stark mixing between S and P states and as was also shown by the pp calculations. The 
results now show slightly more sensitivity to the strong interaction parameters Tis and r 2 p. 
However the fraction of E~ decays largely stays unaltered. 
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Using S and P state capture fractions obtained in this way, with values for Vis and T 2 p 
given in Table II of Ref. | 30|| , the probability for the production of the H particle from 



atoms can be calculated J36|. Aerts and Dover JHJ also give in Fig. 4 of their paper, values, 
as a function of the mass of the H particle m#, for the branching ratios Rs and Rp defined 
by, 

R s = r((S~d) ls - Hn)/T ls (7) 

and 

R P = T((E-d) 2P -> Fn)/r 2P . (8) 

Weighting the fractions of S and P state capture by the corresponding Rs and Rp (values 
for model D in Ref. |3(J were used), then gives the total probability for H particle production 



as a function of the binding energy B(H). Here B(H) is the binding energy of the H particle 



relative to the AA mass. The results are shown in Fig. [12] for four different values of ksrx 
= 0.5, 1.0, 2.0 and 5.0 (curves labelled 1, 2, 3, and 4 respectively) and for K = 1.0 and 7.6 
(dashed and solid lines respectively). 

The results of |3IJ show that H particle production is more likely to proceed from S states 
rather than P states by a factor varying from 1.4 at niu = 2.23 GeV (B(H) = MeV) to 14 
at m H = 2.18 GeV (B(H) = 50 MeV). As a result the probability of H particle production 
increases with Stark mixing (increasing ksrx or K ) as can be clearly seen from Fig. |12[ 

On the other hand the branching fractions Rs and Rp both decrease as the binding energy 
B(H) increases. However the fraction of S state capture increases somewhat from 0.46 to 
0.57 and the P state fraction decreases from 0.38 to 0.27 as B(H) increases from to 50 MeV. 
These changes are largely due to the reduction in the 2P state strong interaction width T 2 p 
from 0.61 to 0.25 eV as B(H) increases; the above capture fractions were calculated for 
ksTK = 1-0 and K = 7.6. As a result, due to the increase in S state capture, the decrease 
in H particle production with increasing B(H) is less rapid than might be deduced solely 
from the overall reduction in the branching fractions for H particle production from the 
atom. 



VI. SUMMARY 

Almost no quantitative information on the interaction of 5~ hyperons with nuclei is 
available at present and it is unlikely that conventional measurements of particle energies to 
investigate H~ hypernuclei will have sufficient accuracy to alter this situation. In contrast, 
the usual precision for measuring the energies of X-rays from transitions between levels of 
hadronic atoms offers the possibility of obtaining quantitative information on the interaction 
of S~ hyperons with nuclei. In the present work we have been guided by the successful 
observation and reasonably precise measurement of strong interaction effects in £~ atoms, 
which has led to our best knowledge, so far, of the interaction of S~ hyperons with nuclei. 

Full atomic cascade calculations have been performed for £~ and S~ atoms and con- 
firmed, as expected, that the processes within these two hadronic atoms are very similar. 
The remaining major differences are in the production reactions. Whereas S~ hyperons are 
produced by the p(K~, 7r + )S~ reaction at rest, the p(K~, K + )S~ reaction occurs at higher 
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energies, thus causing decay losses during the slowing down time of the S~ particle to be 
non- negligible. Prior to such an experiment it will be necessary to optimize the experimental 
setup, which includes a hydrogen production target, a heavy moderator such as Pb or W, 
the target to be studied and the detectors, both for X-rays and for the detection of the 
outgoing K + , which is essential in order to reduce background. 

In the present work we have confined ourselves to studying the dependence of strong 
interaction shifts and widths in S~ atoms on the various parameters of the problem, including 
the nuclear charge. Adopting an optical potential for the strong interaction between the S~ 
and the nucleus and assuming it to be attractive and 15-20 MeV deep with an imaginary 
part of 1-3 MeV, we are able to propose four targets along the periodic table, namely, F, 
CI, I and Pb as suitable for X-ray measurements. This also proves to be a sensible choice of 
targets for moderate changes in V op t. Nevertheless, even if the actual potential turns out to 
be very different from the one used in the present calculations, for example using potential 
G of Sec. [TV], only relatively small changes will result, because the whole phenomenon of 
hadronic atoms is dominated by the Coulomb interaction. Starting from a G-matrix based 
EN interaction, we constructed a version of (potential G) which would indeed lead to a 
slightly different choice of optimal targets. It was shown that the qualitative features of such 
a potential can be deduced by an inversion procedure, where one fits a phenomenological 
DD potential to the strong interaction data deduced from S~ atom experiments. Although 
we do not consider potential G as a realistic one, this example served to demonstrate the 
feasibility of the inversion procedure. 

Finally, the special case of E~d atoms has been discussed in great detail, in view of its 
role in current experiments aimed at the H dibaryon. The role of Stark mixing, its effect 
on S and P state capture of S~ by the deuteron, together with estimates of the resulting 
probability for producing the H dibaryon have been considered in detail. 



Some of the early calculations on H d atoms were carried out in collaboration with G. T. 
A. Squier. The estimates of H particle production from E~d atoms were made by J. Lowe. 
This research was supported by THE ISRAEL SCIENCE FOUNDATION founded by The 
Academy of Sciences and Humanities. 
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FIG. 1. Calculated absolute yields of X-ray transitions in £~ Pb atoms as a function of the 
initial population parameter a. Results are shown for raj — > nj transitions with I = n — 1. 
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FIG. 3. Calculated strong interaction widths and upper level relative yields for the 9K level 
heavy H~ atoms. 
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FIG. 4. Solid curves: calculated strong interaction widths and upper level relative yields for 
the 71 level in medium-heavy H~ atoms. The dashed curves are for bo = —0.25 + i0.04 fm, i.e. a 
repulsive real potential. 



17 



H atoms 4F states 



1.00 



0.80 



| 0.60 
S 



a. o.4o 



0.20 



0.00 



13.0 14.0 15.0 16.0 17.0 18.0 19.0 20.0 21.0 

10.00 r 1 1 1 1 1 1 1 



> 
CD 



1.00 



0.10 



0.01 1 1 1 1 1 1 1 1 

13.0 14.0 15.0 16.0 17.0 18.0 19.0 20.0 21.0 

atomic number 

FIG. 5. Calculated strong interaction widths and upper level relative yields for the 4F level 
atoms. 
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FIG. 6. Calculated strong interaction shifts and widths for the 4F level in E~ atoms of CI for 
different optical potentials, see text for detail. 
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FIG. 7. Calculated strong interaction widths and upper level relative yields for the 3D level 
atoms. 
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FIG. 8. Calculated X-ray spectrum for a teflon target. 
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FIG. 9. Real optical potentials for S interaction with Sn. Solid curve (G) is for Eq. @ used 
to generate the 'pseudo data'. The other curves are the fitted potentials obtained under various 
constraints. See text. 
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FIG. 10. Calculated fractions of S and P state capture, and of H~ decays, for three choices 
of the strong interaction parameters T\s and T2P (see text) as a function of the Stark mixing 
parameter ksTK- 
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FIG. 12. Calculated probabilities for H particle production as a function of its binding energy 
B(H). Curves labelled 1, 2, 3 and 4 correspond to values of the Stark mixing parameter ksTK = 
0.5, 1.0, 2.0 and 5.0 respectively. The dashed and solid lines correspond to values of Kq = 1.0 and 
7.6 respectively. (See text). 
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TABLES 

TABLE I. Calculated binding energies and strong interaction widths (in keV) of 3D E~ atomic 
states. 

™c ii^ ii5 

Coulomb only 126.36, L=0.000 174.71, L=0.000 230.90, L=0.000 

potential 1 126.39, L=0.012 174.81, L=0.052 231.32, L=0.226 

potential 2 126.38, L=0.010 174.78, L=0.040 231.17, L=0.167 



TABLE II. Predictions for likely targets for a S atoms experiment. Calculations are based on 
a tp potential with bo = 0.25 + i0.04 fm. ~E X is the transition energy, Y is the upper level relative 



yield. 

target F CI Sn I Pb 

transition 4F 3D 5G -»■ 4F 8 J -> 71 8 J -> 71 VOL -> 

E x (keV) 131.29 223.55 420.25 474.71 558.47 

y 0.31 0.37 0.76 0.43 0.58 

shift (keV) 1.56 1.84 0.67 2.79 1.73 

width (keV) 0.99 1.14 0.43 2.21 1.26 



TABLE III. Parameters for the G-matrix interaction Eq. (Q) 

Pi (fm) 2^0 09 0.5 

a. (MeV) -2.595 -200.7 821.4 

h (MeV fm) 0.2632 63.94 -124.0 



TABLE IV. Calculated shifts and widths (in keV) due to potential G, Eqs. (g,|) and Table [TfiJ. 



atom 


0(3D) 


F(3D) 


S(4F) 


C1(4F) 


Sn(7I) 


1(71) 


W(9K) 


Pb(9K) 


shift 


3.343 


14.23 


7.511 


20.53 


8.725 


40.39 


3.191 


22.90 


width 


0.924 


7.619 


2.523 


11.68 


3.393 


40.29 


0.846 


16.63 
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